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Abstract 

I review the status of today's theoretical description of Standard Model pro- 
cesses relevant for Tevatron and LHC analyses, and of the tools that are used 
in phenomenological studies. I will also discuss a few recent ideas to further 
refine our abilities to perform technically challenging calculations. 

CN ■ 1 Introduction 

t— i ■ 

Today we face many fundamental questions, some of which are driven by experimental data, such as 
the question about the mechanism of electroweak (EW) symmetry breaking, the nature of dark matter, 
and the physics associated with the vacuum energy, as well as questions that are driven by theoretical 
^ | curiosity and ambition, which are in essence propelled by our hope that there is an elegant structure 

^-j- ■ behind what we observe in nature. Questions of the latter type include why there are three generations, 

what causes the hierarchy of fermion masses and mixing, and how the strong CP problem is resolved. 
While questions of the first type do have a definite answer, this is not necessarily the case for the questions 
1 of the second type. From this point of view, some of these questions might not even be the right ones 

^ to ask. With the start of the LHC we feel that we are at a verge of big changes, the depth of which we 

Q-i can not assess yet. Indeed, the LHC marks the start of a long research program and experiments at the 

LHC are expected to revolutionize our understanding of the fundamental forces and matter. The LHC 
will definitely explore the origin of mass and the associated nature of EW symmetry breaking. In the 
. course of this, it might also shed light on the nature of dark matter and the origin of the matter-antimatter 

£> | asymmetry. It may also explore the physics that underlies the evolution of the early universe. While it is 

clear that the LHC will not answer all the fundamental questions that we have, the questions we ask now 
will most likely change after the LHC era. It is therefore really a great and unique time to be a particle 
\ physicist. 
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Since the beginning of its run in 2010, the LHC has been remarkable successful. ATLAS and CMS 
collected around 45 pb _1 in 2010, more than 1 fb _1 by the time of summer conferences, and more than 
5 fb _1 by the end of the year. At the end of the 201 1 run, almost every week has set a new record in in- 
stantaneous luminosity. With the 2010 and early 2011 data, remarkably, all major Standard Model (SM) 
processes have already been re-established, including single-top and di-boson production, challenging 
measurements (because of the small/large cross-sections/backgrounds) that have been performed at the 
Tevatron only in recent years. By now, we have entered a new territory in the search of physics beyond 
the SM (BSM) with sensitivities already well exceeding those of LEP and the Tevatron. 

One important question then concerns the role of QCD for LHC measurements and new-physics 
searches. Understanding how QCD works is essential in order to make accurate predictions for both 
the signal and background processes. This typically requires complex calculations to higher orders in 
the perturbative expansion of the coupling constant. Understanding QCD dynamics can however also 
help reduce backgrounds and sharpen the structure of the signal. This can for instance be achieved by 
designing better observables, by employing appropriate jet algorithms, by using jet-substructure, or by 
exploiting properties of boosted kinematics. Finally, once discovery is made, QCD will be crucial to 
extract the properties (masses, spins, and couplings) of the new states found. Therefore, at the LHC, no 
matter what physics you do, QCD will be part of your life. 

It is interesting to first recall a recent measurement that was the origin of considerable excitement. 
In April 2011, CDF reported the observation of a peak in the rrijj distribution in W + dijet events ID, 




Fig. 1: The dijet invariant mass distribution for W plus dijet events as measured by CDF. The left pane shows the 
fits for known processes only. The right pane shows, by subtraction, the resonant contribution to rrijj including 
WW and WZ production and a hypothesized narrow Gaussian contribution. Figures taken from [ 1]. 



see FigQ] The first measurement had a 3.2a significance, and was based on 4.3 fb -1 . Subsequently, 
more data (7.3 fb _1 ) has been analyzed, leading to a significance of more than 4cr [2]. Since then, a large 
number of tentative new-physics explanations appeared on the arXiv, along with a few SM analyses 
that address the question of whether this effect can be attributed to a mismodelling of one of the SM 
backgrounds (in particular single top) 0431. The excitement was curbed few months afterwards, when in 
June, DO announced that it did not confirm the excess seen by CDF [6]. It is yet unclear what the reasons 
for the discrepancy between CDF and DO findings are, if any. However, this example demonstrates 
that even in the case where one identifies a mass peak in the tail of a distribution (a scenario that was 
considered "an easy discovery") a robust control of SM backgrounds remains mandatory, in particular 
when the shape of the backgrounds is one of the issues. Currently we have a number of other recent 
measurements at collider experiments that report a few sigma deviations from the SM predictions. This 
is for instance the case for the top forward-backward asymmetry measured by both CDF [7] and DO (SJ, 
for the dimuon charge asymmetry measured by DO [9], for W + b measured by DO ifTOl . CP violation in 
time-integrated D° decays [11], and a few more. 

The important question becomes then what the tools at our disposal are to make precise predic- 
tions, and whether we have the solid control of signal and background processes that is needed in order 
to claim discoveries. In this write-up I will discuss the status of our theoretical knowledge of most 
important SM processes, the tools at our disposal to describe these processes, and the impact of QCD 
higher orders in view of recent Tevatron/LHC results. Finally, I will discuss a few recent ideas to further 
improve on the way we perform technically challenging calculations. 

2 Perturbative tools 

The range of physics analyses that one can do at the LHC is very broad. It includes pure instrumental 
QCD studies, such as measurements of partem densities and inclusive jet cross-section measurements, 
precision EW measurements, Higgs searches, direct and indirect BSM searches, B physics, top physics, 
diffractive studies and forward physics, and heavy ion physics. Each of these topics includes a vast 
number of measurements and studies. Yet, there are three things that everybody involved in any of these 
analyses can not live without. These are Monte Carlos (MCs), partem distribution functions (PDFs), and 
jets. 
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Fig. 2: The multiplicity of jets with pt > 55 GeV for events in a W enriched control region at the LHC (7 TeV). 
Figure taken from lfT8ll . 

2.1 Monte Carlos and leading order matrix elements 

The first thing one can not live without at the LHC are MC generators. Apart from very few excep- 
tions, every analysis at the LHC uses a MC program for the simulation of the signal process, for the 
backgrounds, for subtracting the underlying event and the non-perturbative contributions, and/or for ef- 
ficiency studies and modeling of the detector response. The current level of sophistication is such that 
essentially not a single study relies on Pythia/Herwig alone. It is well understood that in multi-parton 
processes it is important to describe the multiple hard QCD radiation at least using exact matrix elements, 
employing for instance Alpgen [12], Madgraph Ifl3l . or Sherpa lfl4lL 

Since experimental studies rely heavily on all these leading-order (LO) tools, there is continuous 
progress in their development, and the Herwig/Pythia codes that we have today bear little resemblance 
to their original version of the '80s. In particular, in Pythia 8.1 lfT31 (a C++ code) there is a new 
fully interleaved p t -ordered multi-parton interaction (MPI), initial- and final state evolution (the original 
mass-ordered evolution is not supported any longer), a richer mix of underlying event processes (7, 
J/vP, DY), the possibility to select two hard interactions in the same event, an x-dependent proton size 
in the MPI framework, the full hadron-hadron machinery for diffractive systems, several new processes 
in and beyond the SM, and various other new features. Herwig++ lfl6l (the current version is 2.5.2) 
has new next-to-leading order (NLO) matrix elements, including weak boson pair production, a colour 
reconnection model, diffractive processes, additional models of BSM physics, and new LO elements for 
hadron-hadron, lepton-lepton collisions, and photon-initiated processes. Sherpa lfT4ll (version 1.3.1) has 
improved integration routines in Comix, a simplified kinematics reconstruction algorithm of the parton 
shower (PS), leading to numerically more stable simulations, HepMC output for NLO events and various 
other improvements/bug-fixes. Madgraph lfL3l (version 5) has a completely new diagram generation 
algorithm, which makes optimal use of model-independent information, has an efficient decay-chain 
package, and a new library for the colour calculations. Some applications from ALOHA [ 17 ], an automatic 
library of helicity amplitudes for Feynamn diagram computations are also implemented. Altogether, 
there is continuous, fast progress in various directions. So far, it is amazing how well these tools work, 
once the normalization is fixed using data. A very recent comparison of data with Alpgen up to six jets 
(a control region for BSM searches) is shown as an illustration in Fig. |2] [ 18]. Yet, the devil is often in 
the detail (i.e. in the ~ 20 — 30% effects that are opposite to expectations). For instance in general one 
expects matrix-element based MCs to work better than pure PSs, but this is not always the case (see e.g. 
Ifl9l ). Altogether, these LO programs will undergo a real stress test in the coming years. 
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2.2 The NLO revolution 

Theorists like to advertise NLO computations by using the reduction of scale uncertainties in the predic- 
tions as an argument, which is meant to reflect the reduction in the theoretical perturbative uncertainty. 
However, the strongest argument in support of NLO calculations is their past success in accurately de- 
scribing LEP and Tevatron data. Because of the importance of NLO corrections, an industrial effort has 
been devoted in the last years to these computations [20]. Recent revolutionary ideas in the way NLO 
computations are performed include sewing together tree-level amplitudes to compute loop amplitudes 
(using on-shell intermediate states, cuts, unitarity ideas, . . . ) [21], the OPP algorithm, an algebraic way to 
extract coefficients of master integrals by evaluating the amplitudes at specific values of the loop momen- 
tum ll22l . and D-dimensional unitarity, a practical numerical tool to evaluate full amplitudes, including 
the rational part, with unitarity ideas [23]. For a pedagogical review on unitarity methods see 1241 . These 
methods led in the past 2 to 3 years to a number of 2 — > 4 calculations at hadron colliders. These include 
W + 3 jetsjp[26l, Z + 3 jets E3, tibb EHJ, ti -> W + W~bb [29], W+W+ + 2 jets (301, W+W~ + 2 
jets Oil , bbbb 021 . tt + 2 jets OBI , four jet production Q4l . and a few other ones. 

Feynman diagram methods have also been applied successfully to 2 — > 4 calculations, this is for 
instance the case for quark-induced bbbb [35], tibb P6l . W + W~bb 071 production, W^^+]et [38], and 
a number of vector boson fusioon (VBF) processes which are available in the public code VBFNLO P9l . 
VBF processes have never been measured at the Tevatron. First measurements of VBF W or Z produc- 
tion at the LHC are therefore of particular interest. They also pave the way to measurements of VBF 
Higgs production. Note that only a few years ago, performing this type of calculation with Feynman 
diagrams was considered an impossible task. 

A novel approach to NLO calculations, based on the results of [40], promotes traditional tree 
algorithms to generators of loop-momentum polynomials, that are called open loops [41 ]. The excellent 
performance of the method has been demonstrated, open loops have therefore a potential to address a 
number of multi-particle processes at hadron colliders. Another recent approach evaluates one-loop QCD 
amplitudes purely numerically Il42l . The algorithm consists of subtraction terms, approximating the soft, 
collinear and ultraviolet divergences of one-loop amplitudes and a method to deform the integration 
contour for the loop integration into the complex space. The algorithm is formulated at the amplitude 
level and does not rely on Feynman graphs. The power of this method has been demonstrated recently 
with the leading color NLO calculation of up to seven jets in e + e~ collisions ll43l . 

Given that both Feynman diagram and unitarity based methods allowed us to compute 2 — > 4 
processes at NLO in QCD, it might be unclear where the revolution advocated in the heading of the 
subsection lies in. The revolution, I believe, is not yet in the applications that we see today, rather in the 
prospect for low-cost fully computer-automated NLO calculations even beyond 2 — > 4 in the near future. 
Indeed, two 2—^5 processes have already been computed at NLO, namely W + 4 jets [44] and Z + 4 
jets (45]Q Fig.[3]illustrates in the case of W~ +4 jets the typical effect of including NLO corrections: one 
obtains a considerable reduction of the scale uncertainty, and, for some distributions, a change in shape. 
As far as the full automation is concerned, let me highlight three interesting approaches. The first one 11461 
is a method based on Feynman diagrams, it uses the OPP procedure for the virtual calculation, and the 
FKS subtraction of divergences, together with clever and efficient procedures to deal with instabilities. 
More improvements and refinements are to be expected soon. At present there is no public code, instead 
the idea is to provide iV-tuples. The second approach, HELAC-1L00P [47], is a program that evaluates 
numerically QCD virtual corrections to scattering amplitudes using the OPP method. The public program 
is part of the HELAC-NLO framework that allows for a complete evaluation of QCD NLO corrections. 
GoSam [48 ] is a third approach which aims at the full automated calculation of NLO corrections for multi- 
particle processes. The one-loop amplitudes are generated using Feynman diagrams and are reduced 
using D-dimensional unitarity, a refined Passarino-Veltman style tensor reduction, or a combination of 

'in both cases the leading colour approximation has been used. This approximation is expected to give rise to very small 
(percent) corrections only. 
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Fig. 3: The transverse momentum distribution of the leading four jets in W~ +4-jet production at the LHC (7 TeV) 
at LO and NLO. The lower panels show the LO and NLO scale-dependence bands (H' /4 < fiR = [ip < Hip) 
normalized to the central NLO prediction (/i^ = /ip = H' T /2). Figure taken from [44 1. 



both. GoSam can be used to calculate one-loop corrections to both QCD and electroweak theory, and 
model files for BSM theories can be also linked. An interface to programs calculating real radiation is 
included too. The flexibility of GoSam has been demonstrated explicitly by considering various examples. 

2.3 Merging NLO and Parton Showers 

While NLO predictions provide relatively accurate results for inclusive cross-sections, they do not furnish 
an exclusive description of the final state that can be compared with actual particles in the detectors, as 
MC programs do. It is therefore useful to combine the best features of both approaches. Two public 
frameworks exist for this purpose, namely MC@NL0 [49] and P0WHEG [50]. These tools are almost 10 years 
old now, and since their conception a long list of processes has been implemented in both frameworks. In 
particular, recently the P0WHEG BOX was released f5T1 . which is a general framework for implementing 
NLO calculations in shower MC programs according to the P0WHEG method. The user only needs to 
provide a simple set of routines (Born, colour-correlated Born, virtual, real, and phase space) that are 
part of any NLO calculation. 

Similarly, aMC@NL0 is a novel approach to a complete event generation at NLO. aMCONLQ has been 
used recently for the calculation of W/Zbb [52] and W+dijet production [53]. Fig. [4] (left) shows an 
application to Higgs searches of the W/Zbb calculation: the invariant mass of the pair of the two leading 
6-jets, for the processes Wbb, Zbb, WH, and ZH. The figure illustrates a case where signals and 
irreducible backgrounds are computed with the same accuracy. The process Wbb has been implemented 
shortly before also in the P0WHEG BOX [54]. 

Fig. H] (right) illustrates predictions from aMC@NL0 for the invariant mass for the dijet system in 
Wjj (the observable mentioned in the introduction where CDF observed a large deviation from the 
SM). CDF and DO estimate the Wjj using a leading order Monte Carlo (LO+PS) re-weighted to the 
NLO cross-section or to data. With aMC@NL0 instead it is possible to compute directly the Wjj cross- 
section at the NLO+PS level. It was therefore particularly interesting to check whether there is any shape 
difference between LO+PS and NLO+PS in the Mjj distribution. The study of ref. ll53l shows that there 
is no sizable shape difference. Another interesting application of aMCNLQ is the calculation of scalar and 
pseudo-scalar Higgs production in association with a ti pair ll55l . 
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Fig. 4: Left: Invariant mass of the pair of the two leading 6-jets for Wbb, Zbb, WH{-^r ivbb), ZH{-+ £ + £~bb) at 
the LHC (7 TeV), the latter two are rescaled by a factor of ten. Figure taken from [52]. Right: Invariant mass of 
the pair of the two hardest jets in W+ jet events with CDF/DO inclusive cuts. Figure taken from ref. Il53ll . 
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Fig. 5: Comparison of NLO and POWHEG+PYTHIA results for the #t,tot distribution in the process W + W + + 2 
jets at the LHC (7 TeV), when all jets are included in the definition of -£/t,tot (left pane), and when only the three 
hardest jets are included (right pane). Figure taken from ref. [56]. 



A lot of effort has been devoted recently also to the implementation of higher multiplicity pro- 
cesses in the PQWHEG BOX and in aMCONLO. The first 2 — > 4 process that has been implemented in the 
PQWHEG BOX is pp ->■ 2 jets including both the QCD induced part (56l as well as the VBF 

contributions ll57l . This is a relatively simple 2 — > 4 process since the cross-section is finite without 
any cut on the jets. As expected, for inclusive observables there are only minor differences between 
pure NLO and P0WHEG+PS, but for exclusive observables, depending on the details of the observable 
definition, there can be important differences. This is shown in Fig. [5] for two different definitions of 
#t,tot = YljPtj' tne transverse energy of the event. From the figure it is clear that if only the three 
hardest jets are included in the definition of Ht, the corrections from the PS are very moderate (right 
pane), but if all soft jets present in the event are included, then additional radiation from the PS can 
alter the distribution substantially (left pane). The P0WHEG method has been used recently also for the 
description of other important processes, e.g. for the calculation of Z boson in association with a top 
anti-top pair at NLO accuracy including PS effects Il58ll . 

2.4 MENL0PS and LoopSim 

In NLO+PS approaches only contributions with one additional jet, relative to the Born contribution, are 
computed accurately, while all other emissions are described only in the shower approximation. There 
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Fig. 6: The LO and NLO distributions for three observables in Z+jet production that are identical at LO: the Z 
transverse momentum (left), the p t of the hardest jet (middle), and the scalar sum of the transverse momenta of all 
the jets, iJxjets (right). The bands correspond to the uncertainty from a simultaneous variation of /i = fin = fip 
by a factor of two around the default fi = Vjjf^ + m 2 z . Figure taken from ll6T1l . 

are however situations in which one has important contributions from higher multiplicity final states 
(e.g. because of new channels that open up that are enhanced by kinematical factors, gluon PDFs etc.). 
MENL0PS [ 59 , 60] is a method to further improve on NLO+PS predictions with matrix elements involving 
more partons in the final state. For example, for W production it includes, as in MC@NL0 or PQWHEG, W 
production at NLO, the PS, but also W + 1, 2, 3, . . . jets using exact matrix elements. Roughly speaking, 
it uses a jet-algorithm to define two different regimes, and then corrects the 1-jet fraction using exact 
matrix elements and the 2-jet fraction using the NLO A-factor. This achieves NLO quality accuracy for 
inclusive quantities but an improved sensitivity to hard radiation and multi-parton kinematic features. 

A further recent theoretical development is LoopSim. If one considers the process W + 1 jet, the 
three observables p t) z, ptj» an d ^Tjets = Y^jPtj w& identical at LO. However, as illustrated in Fig. [6j 
at NLO pt,z has a moderate A-factor (;$ 2), p t ,j has a large if -factor (~ 5) and -ffrjets has a giant 
A'-factor (~ 50). The very large A-factors in the last two observables are due to the fact that the NLO 
result is dominated by configurations where there are two hard jets and a soft W (these are enhanced 
by EW logarithms), additionally there is an important enhancement coming from incoming qq channels. 
LoopSim [61] is a procedure that uses a sequential algorithm, close to the Cambridge/Aachen one, to 
determine the branching history, "loops" over soft particles (i.e. they are removed from the event and 
the residual event is adjusted), and it uses a unitary operator to cancel divergences. In essence, this is a 
way to extend a calculation that is exact at a given order in perturbation theory, in an approximate way 
to higher orders. The procedure is expected to be more accurate the larger the corresponding A'-factor 
is. One might expect other extensions of the MLM/CKKW matching procedure along the same lines as 
MENLDPS and LoopSim in the near future. 

3 Top-quark production 

The top is the most interesting SM quark. Its large mass implies a large Yukawa coupling, which causes 
the top to be a prominent decay product in many BSM models. LHC data have already been successfully 
compared to approximate next-to-next-to-leading order (NNLO) predictions Il62ll63l . however various 
approximate NNLO predictions, based on a threshold resummation, do not fully agree within quoted 
uncertainties Il64l468l . Significant improvements in the tt cross-sections can be expected only upon in- 
clusion of the complete NNLO corrections. A better perturbative control of the top-quark pair production 
cross-section is also important to further constrain gluon PDFs, to have an accurate extraction of the top 
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mass from the cross-section, and to improve our perturbative control over the ti forward-backward asym- 
metry. In fact, an almost 3a deviation from the SM is observed by CDF, which becomes a 4.2a effect in 
the high-mass region, M t i > 450 GeV 0. The large inclusive asymmetry has been seen both by CDF 
and DO [ 8 ], while the rise in the spectrum of the asymmetry is not confirmed by DO. One has however to 
bear in mind that ti production is a difficult measurement given the presence of neutrinos in the final state, 
the combinatorics in the reconstruction of the tops, and the limited statistics at the Tevatron. Neverthe- 
less, various suggestions have been made recently to explain the asymmetry in terms of BSM physics, but 
all proposals face the problem that they have to preserve the good agreement with the symmetric ti ob- 
servables, respect dijet bounds and/or must evade the stringent limits on like-sign top production. Fervid 
activity is therefore currently devoted towards a complete NNLO calculation of ti production (see Il69l 
and references therein). Recently, the program Top++ was released that evaluates numerically the total 
inclusive cross-section for producing top quark pairs at hadron colliders. It calculates the cross-section 
at fixed order through approximate NNLO and it includes a soft-gluon resummation in Mellin space at 
next-to-next-to-leading logarithmic accuracy [70]. The program Hathor irTTl also calculates the total 
cross-section for top-quark pair production in hadronic collisions. It includes approximate next-to-next- 
to-leading order perturbative QCD corrections. It also offers the possibility to obtain the cross-section 
as a function of the running top-quark mass. A direct comparison between the two programs is however 
not possible as they use a different subleading terms in the threshold expansion. 

Besides the inclusive cross-section, top pair production in association with light jets, heavy par- 
tons, photons or heavy vector bosons is also interesting. We know now at NLO ti + 2 jets 11331 . 
tibb E8ll36l tiZ EH, ti + 7 [73], and for ti + 1 jet including hard jet radiation by the top quark 
decay products 11741 . A matching to parton shower has been recently achieved also to ti + 1 jet produc- 
tion 117511761 . 

The NLO calculation of single top in the s and t channel H77j|7H and in the Wt channel 11791,18011 
have also been matched to parton showers. For a discussion of recent theoretical progress in single-top 
physics at hadron colliders, and of aspects of single-top production in beyond Standard Model scenarios 
I refer the reader to ref. [81]. 

4 A few theoretical issues in Higgs production 

Possibly the most awaited result during the first year of running of the LHC concerned Higgs searches. 
Currently, both ATLAS and CMS reached the expected sensitivity, around or better than the SM cross- 
section. ATLAS restricted the most likely mass range at 95% confidence level (CL) to the region 115.5 — 
131 GeV. They observe an excess around 126 GeV, with a local significance of 3.6a, with contributions 
from the main channels H — > tt, H — > ZZ ->• 41, H — > WW — > 2l2v. The global significance 
taking into account the look-elsewhere-effect is 2.3a. CMS excludes the region 127 — 600 GeV at 95% 
CL (while their expected exclusion is 117 — 543 GeV). They could not exclude the region below 127 
GeV since data has a modest excess of events between 115 and 127 GeV that appears, quite consistently, 
in five independent channels. The excess is compatible with a SM Higgs hypothesis in the vicinity of 
124 GeV or slightly below, but the statistical significance, 2.6a local and 1.9a global after correcting for 
the look-elsewhere-effect in the low mass region, is not large enough to say anything conclusive. For 
both experiments, what is observed now is consistent either with a background fluctuation or with a SM 
Higgs boson. More refined analyses and additional data in 2012 will definitely give an answer. 

At the LHC, the Higgs is mainly produced via an intermediate top loop in gluon-gluon fusion. 
The urge to understand the EW symmetry breaking led in the past years to the computation of the most 
advanced theoretical predictions for this process. For instance we now know the main gg — ^ H produc- 
tion mechanism including NLO corrections with exact top and bottom quarks in the loop 1821 . NNLO 
corrections in the large mt limit |[83l485ll . EW corrections [86], mixed QCD-EW corrections [87], and 
resummation of large logarithms possibly with N 3 LO soft effects [88-93 ]. Threshold corrections to the 
boson rapidity distribution are also known 11941 . these higher order corrections stabilise the theoretical 
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predictions under scale variations. Furthermore, the most advanced codes |[95ll96l allow for fully ex- 
clusive decays of the Higgs to 77, W + W~ — > e + ve~v, and ZZ — > 4/. A similar accuracy has been 
reached recently also in associated VH production where, because this process is an important one if 
the Higgs is light, the decay of the Higgs into bb has been considered 1971 . NLO EW corrections to 
WH/ZH production including the vector-boson decays are also known [98 ]. As expected, the EW cor- 
rections, which are at the level of (5-10)% for total cross sections, increase with increasing transverse 
momentum cuts. For instance, for pr,H > 200 GeV, which is an interesting range at the LHC, the EW 
corrections to WH production are of the order of -15% for Mh = 120 GeV. 

The fully differential decay of a light Higgs boson to bottom quarks at NNLO in perturbative QCD 
has been computed recently in |[99l . From a technical point of view, it is interesting to note that this work 
constitutes the first physical application of a novel method of non-linear mappings for the treatment of 
singularities in the radiative processes which contribute to the decay width. The program iHixs II 1021 
has also been recently released, it computes the inclusive Higgs boson cross-section, including QCD cor- 
rections through NNLO, EW collections, mixed QCD-electroweak corrections 1571 . quark-mass effects 
through NLO in QCD, and finite width effects for the Higgs boson and the heavy quarks. Furthermore, it 
allows the evaluation of the cross-section in modified Higgs boson sectors with anomalous Yukawa and 
EW interactions as appearing in some extensions of the SM II1011 I1021. 

Given the high accuracy with which gluon-gluon-fusion has been computed, it is interesting to ask 
what is the actual theoretical uncertainty on this process. Unfortunately, there is today no full consensus 
on this question. Some more conservative estimates quote errors of the order of 40% 111031 (at the Teva- 
tron) and similar uncertainties at the LHC, while other studies suggest that the perturbative uncertainty is 
considerably smaller. Assigning a correct theoretical error is very important when claiming an exclusion 
or an excess, and, at a later stage, when making measurements of the Higgs-boson couplings, which is 
the only way to identify the precise nature of the Higgs boson and EW symmetry breaking. Yet, even for 
the main Higgs-production channel there are still some controversies and subtleties. Most controversies 
have to do with how different sources of errors should be combined, others concern the question of how 
to assign/interpret the perturbative uncertainties. I will illustrate here just two of these issues. 

The soft logarithms appearing in cross-sections can be resummed using an effective theory ap- 
proach. Performing such a calculation requires an introduction of a matching scale, where the full and 
effective theory amplitudes must agree. It is well-known that choosing a time-like (i.e. complex) match- 
ing scale effectively resums it 2 enhanced terms. In 111041 it is suggested that this procedure improves the 
convergence of the perturbative expansion significantly, and reduces the uncertainty of the perturbative 
(NNLO) prediction. This approach is criticized in 11 10511 with the arguments that ir 2 are just numbers, 
so that there is no formal limit in which they dominate, and that only one class of it 2 terms is resummed 
(those that arise from the gluon form factor), but not all of them. In this context, one has to mention that 
perturbative QCD is often about pushing approximations beyond their formal limit of validity, and that a 
given approach should be judged by seeing how well it fares in practice. 

The second issue I would like to mention here has to do with a jet-veto in Higgs searches. As 
can be seen from Fig. |7J in Higgs searches one needs to impose a jet-veto to suppress the large top 
background. Higgs production is then studied in 0-, 1-, 2-jet bins separately in order to maximize the 
sensitivity. Currently, ATLAS uses £>t,veto = 25 GeV, while CMS employs pt,veto = 30 GeV in their 
Higgs searches. 

In H 1 0711 the inclusive, NNLO Higgs production cross-section at the Tevatron is split into 0-,1-jet 
exclusive, and 2-jet inclusive components 

^21 = 66.5%^g(0-jet) + 28.6%^g(i-jet)+4.9%^g(> 2-jets) = [-14.3% + 14.0%] . (1) 

<7tot 

The errors denote the scale uncertainty that is obtained by varying the renormalization and the factoriza- 
tion scale together around a central value mn = 160 GeV by a factor of two. In an NNLO calculation of 
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the data and the background expectation from MC, with the yellow band indicating the total systematic uncertainty 
in the normalization (but not the shape) of the various components. The signal is shown for m# = 150 GeV. 
Figure taken from 1 106]. 



inclusive Higgs production, only the 0-jet bin is known at NNLO, while the 1-jet bin is know at NLO and 
the 2-jet bin is computed at LO only. Therefore is it not surprising that the relative errors increase with 
the number of jets. In 111071 one can also find a detailed discussion of why it is not appropriate to use 
the standard scale variation as an estimate of the perturbative uncertainty for the 0-jet bin cross-section 
(Fig.l of [108 ] also shows that for a particular choice of pr,veto one obtains a vanishing scale uncertainty 
band in the 0-jet bin). The numbers in eq. CD were updated by Campbell et al. in [109] who evaluated 
the 2-jet bin contribution at NLO. The effect of this addition was a slight change in all relative numbers, 
and, mainly, a decrease in the perturbative uncertainty of the 2-jet bin, 

— = 60%lg(0 - jet) + 29%+2g(l - jet) + ll%+^%(> 2jets) = [-15.5% + 13.8%] . (2) 

<7tot 

From eq. ©, it is evident that the scale uncertainty is smaller for the exclusive measurement with 0-jets, 
than the one of the fully inclusive measurement. To explain this feature, Stewart and Tackmann recall 
that there are two mechanisms at work in the 0-jet cross-section II110II : there is a large isT-factor from 
perturbative higher orders, as well as large negative logarithms — a s CA/vr In 2 Mn/pt,veto that become 
more important the smaller pt.veto is. They therefore suggest that the error on the 0-jet bin should be 
computed taking into account a full correlation between jet-bins, i.e. the error from the 0-jet cross- 
section is computed from the relation un = 0"i nc i — <7>i_j c t- One obtains then simply A 2 <7n = A 2 <7 inc i + 
A 2 <7>i_j et . The effect of this is illustrated in Fig. [8] While this procedure is certainly more conservative 
than a conventional scale variation, it is clear that to reduce the uncertainty on the jet- veto cross-section, a 
resummation of large logarithms involving the ratio pt,vcto/MH is required. Currently, only resummation 
for quantities related to the jet- veto exist, e.g. for PT.ffiggs 111 1 11 or for the beam-thrust 111 121 . Both 
observables are however not the ones used in current Higgs searches. Furthermore the beam thrust has 
the drawback that it receives very large non-perturbative corrections, as can be easily seen by running a 
PS program once at parton and once at hadron level. 
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5 Gauge boson production processes 

The physics programme involving gauge bosons is particularly rich. The road towards precision mea- 
surements and searches starts from the measurement of inclusive W and Z cross-sections, which, from 
a theoretical point of view, are the most precisely known processes at hadron colliders. Beyond mea- 
surements of purely inclusive W/Z production cross-sections, it is possible to study the production of 
W or Z bosons in association with one or more jets, and ratio of cross-sections. Interesting ratios are 
for instance a(V + (n + 1) }ets)/a(V + njets), with V = W, Z, that start at 0{a s ) in the perturbative 
expansion in the coupling constant, or ratios of a(W ± + njets) /a(Z + njets), that are of order 0{a Q s ). 
In both cases, one expects many experimental and theoretical uncertainties (e.g. those related to the 
choice of renormalization or factorization scale, or uncertainties in the parton disuibution functions) to 
largely cancel in the ratio. A further extension of simple ratios are asymmetry distributions, for instance 
the W or lepton charge asymmetry. These distributions provide strong constraints on parton distribution 
functions and are useful probes of new physics. Other interesting observables measure gauge bosons 
produced in association with heavy quarks (charm or bottom quarks). These processes are particularly 
interesting because of the discrepancies between theoretical predictions and Tevatron data [10], however 
the perturbative calculation of these processes and of the related theoretical uncertainty is challenging. 
Finally, di-boson cross-sections are sensitive to any type of new physics that would modify the trilinear 
gauge couplings, and would give rise to so-called anomalous gauge couplings. These measurements 
are complementary to the ongoing direct searches for beyond the SM physics, and are able to probe 
new-physics scales that are not directly accessible. 

Fig. [9] gives the cross-section for main processes involving gauge bosons at the LHC. The figure 
illustrates that with lfb 1 ATLAS and CMS could collect 0(1O 6 ) and O(10 5 ) W and Z events per 
experiment and per lepton channel. One also sees that, including all lepton channels, lfb" 1 of data 
contain about 100 WW and 10 ZZ events. This means that even with the data available after a first year 
of running, a number of interesting analyses could be performed. 

5.1 Drell-Yan 

The most important gauge boson production process is Drell-Yan. This is the best known process at 
the LHC: it has been computed through NNLO in QCD, fully differential in lepton momenta including 
spin-correlations, EW corrections, finite- width effects, and 7* jZ interference. State-of the art codes are 
described in 11114111151 . Calculations to all-orders also exist, for instance the next-to-next-to-leading log- 
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Fig. 9: Production cross-sections for processes involving gauge bosons at the LHC (at y/s — 7 TeV). Figure taken 
from 01 131. 



arithms (NNLL) transverse momentum resummation [116] and soft gluon resummation have been com- 
puted 11171 . Threshold enhanced corrections to rapidity distributions of Z and bosons at hadron 
colliders are also known at N 3 LO level [118]. These accurate perturbative calculations have been avail- 
able for some time, and now that precise LHC data has been compared to those predictions, one can 
not but praise the impressive agreement between NNLO theory and experiment (see e.g. Fig. [10] and 
ref. 11191 ). One thing to note is that in these comparisons of theory with experiment, the dominant error 
is the theoretical one, however this is mainly dominated by the luminosity uncertainty (of the order of 
4%). 

At this level of precision, it is legitimate to start worrying about mixed QCD and EW correc- 
tions. QCD+EW interference is formally a higher-order effect, and exact results are not known currently. 
However the dominant effects come from emissions in the soft/collinear regions where an approximate 
factorization holds. It turns out that in general mixed corrections are small, O(0.5%), but they are en- 
hanced in the region of large transverse momenta, where they can reach up to around 5% 11211 . 



5.2 Charge asymmetry 

The natural extension of the inclusive cross-section is the Ryy = W + jW~ ratio. One can then study 
Rw as a function of kinematical variables, e.g. one can look at the charge asymmetry as a function of 
lepton rapidity rj 

Rwiv) 



1 



(3) 



Rwiv) + 1 

This measurement is very sensitive to PDFs since in the ratio asymmetric properties of PDFs are en- 
hanced, while many uncertainties cancel. Fig. [TT] shows the relative good agreement of theoretical pre- 
dictions that use various PDFs with ATLAS data. It also illustrates how the shape of the theoretical 
prediction is sensitive to the PDFs chosen. Indeed, ATLAS and CMS measurements of this distribution 
have been already used by the Neural Network (NN) collaboration to constraint PDFs. In particular a 
reduction of uncertainty of the order of 10 — 30% in the range x = [10~ 3 , 10 _1 ] was obtained for the 
valence- and sea-quark distributions. It is interesting to observe that LHCb data at larger rapidities probe 
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Fig. 10: Comparison of NNLO theory and CMS data for Drell-Yan observables. Figure taken from ref. II 1 2011 . 
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Fig. 11: Predictions for the W lepton asymmetry at NLO, obtained with DYNNLO 01141 using the CT10 111221 . 
MSTW08 1 123 1 and NNPDF2.1 II 1241 parton sets, compared to measurements for the muon charge asymmetry 
from ATLAS Ifl25l (7 TeV). Figure taken from ref. |[T26l . 



larger and smaller values of x that are currently less constraint. They will therefore soon have a larger 
impact in PDF determination than ATLAS and CMS have. 

5.3 Production of W/ Z boson in association with jets 

At the LHC, because of the large energy, the production of W and Z bosons in association with jets is 
very likely. This is illustrated in Fig. [12] which shows the differential distribution for H^, the total trans- 
verse energy of the event, for various jet multiplicities. Since the cross-sections with an additional jet is 
rescaled in the figure by a factor 1CP 1 compared to the cross-section with one less jet, it is evident that at 
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Fig. 12: Cross-section for W — > ev and jets at the LHC (7 TeV) as a function of the total transverse momentum of 
the event H^. Figure taken from II 1271 . 

high Ht; (a region of particular interest for various new -physics searches), all jet-multiplicities contribute 
similar amountsll Because of this, it becomes very important to have a good perturbative control of pro- 
cesses involving the production of W/Z bosons together with many jets. The perturbative calculation of 
processes involving a large number of jets is quite difficult beyond LO. However, as discussed in Sec. 12.21 
recent years have seen a revolution in the techniques used for NLO calculations. These novel techniques 
allowed in the last five years the calculation of a large number of processes involving gauge bosons and 
jets. In particular, while V + l and V + 2 jets have been described to NLO in QCD since 1983 II128II1291 
and 2002 lH5TM3"2l . to quote a few examples we know now at NLO VV + 1 jet ifBBTIBol W + 3 
jets (25l|26l, Z + 3 jets E2, W+W+ plus dijets ffl , W + W~ plus dijets [31 ] W + W~bb [|29j|37] and 
W + 4 jets Il44l and Z + 4 jets P31 . Furthermore, various vector boson fusion (VBF) induced gauge 
boson production processes have been computed and are available in the public code VBFNLO [39]. 

Fig. [13] shows the transverse momentum distribution for the four p t ordered jets in Z + 4 jets 
production. The middle pane illustrates the typical reduction of the dependence of the cross-section on 
renormalization and factorization scale at NLO, compared to LO. It is however also evident that while for 
some distributions (e.g. Pt,j4) the ratio of LO/NLO is flat, in other cases the shape of the NLO distribution 
is different from the LO one, so that LO/NLO has a non-flat slope. Is it also interesting to look at the 
ratio of Z/W + and Z/W~ displayed in the lower pane. The agreement of LO and NLO predictions 
for these ratios illustrates the excellent perturbative control that one can achieve at NLO on these ratios. 

2 Of course, this statement depends on the precise definition of the jets, in particular on their p t cut. 
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Fig. 13: Transverse momentum distributions of the leading four jets in Z + 4-jet production at the LHC (7 TeV). 
Figure taken from l45ll . 

Furthermore, the fact the Z/W + is flat, while Z/W~ rises with the p t can be attributed to the fact that 
both W + and Z production are dominated by the u-quark distribution, while W~ is mainly produced 
from d quarks in the protons, therefore the Z/W~ ratio displays the slope of the u/<i-quark distributions 
as a function of the jet transverse momenta. 

5.4 Di-boson production and anomalous couplings 

The ZZ, WW, WZ di-boson production of processes that have been implemented recently in the 
P0WHEG BOX. The calculation includes ^y*/Z interference, single resonant contributions, interference 
effects for identical fermions and, for WW and WZ the effect of anomalous couplings. The gluon- 
gluon fusion contribution, that is formally NNLO, but is important when Higgs search cuts are applied, 
is available the gg2WW and gg2ZZ generators 111371 - 11391 , and, recently, also from the program MCFM H1401 . 

A pure NLO calculation, as implemented earlier in MCFM II 14 1 II 14211 . reveals that for these processes 
the conventional scale variation of the LO result is very modest but underestimates completely the size 
of the NLO corrections. E.g. for pp — > W + W" — >• e + u e fi~u ll production at the 7 TeV LHC without 
any cuts, the LO cross-section using NNPDf2.1 111241 is 375.2lg| fb, while the NLO cross-section is 
499.81$ fb Ifl43l . Here the error denotes the scale uncertainty. Similar results (with smaller cross- 
sections) are obtained for the other processes. The reason for the large NLO corrections (not caught 
by the LO scale variation) is that new partonic channels open up at NLO. It is therefore clear that only 
a NLO calculation can provide a reliable estimate of the cross-section and of its error. These di-boson 
production processes are particularly interesting since they are important backgrounds to Higgs searches. 
Furthermore, they are sensitive to new physics at high scales through the measurement of anomalous 
trilinear couplings (ATGCs). Indeed, while the LHC does probe new physics at the TeV scale directly, 
ATGCs indirectly probe physics in the multi-TeV range, since they arise when high-energy degrees of 
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Fig. 14: The integrated cross-section for W + W~ production at the LHC (7 TeV) as a function of the cut on the 
transverse momentum of the leading lepton is shown on the right. Figure taken from [ 14-3 1 . 



freedom are integrated out. Both the Tevatron 11144141461 and LEP 111471 were able to place quite stringent 
bounds on ATGCs. However, since their effects are enhanced at high energies, one expects even better 
bounds from the LHC. Indeed, CMS already presented bounds on the anomalous couplings appearing in 
an effective Lagrangian with the parametrization of ref. 111481 without form factors 111491 . 

Following refs. 11148 [I1501I151H . one can parametrize the most general terms for the WWV vertex 
(V = 7, Z) in a Lagrangian that conserves C and P as 



Cos = igwwv gX(W*w»v v - w^w^v") 



■Vi.,.. , 

\V (4) 

+ k v W*W v V» v + —W*„WnV pl " , 



m, ^ p 



where W^ u = d^W u - dyW^, g wwz = -ecot6 w and g WWl = -e. In the SM g\ = k\ = 1, and 
A^ = 0. Any departure from these values (AgY = g\ — 1 etc.) would be a sign of new physics. In the 
PQWHEG generator, all six parameters can be set independently. 

In the presence of anomalous couplings, the effective Lagrangian of eq. ((U) gives rise to interac- 
tions that violate unitarity at high energy. Thus, in order to achieve a more realistic parametrization, the 
couplings are multiplied by form factors, that embody the effects arising from integrating out the new 
degrees of freedom. The precise details of the form factors therefore depend on the particular model 
considered. Paralleling the discussion of ref. 111521 . in the P0WHEG BOX it is assumed that all anomalous 
coupling Ag are modified as 

(1 + M2 V /A2)2' (5) 

where M vv is the invariant mass of the vector boson pair and A is the scale of new physics. 

Fig. [14] sets the anomalous coupling to the maximum deviation from the SM allowed by LEP 
bounds and displays the sensitivity of the transverse momentum of the leading jet in W + W~ events to 
anomalous couplings, for the case of a form factor A = 2 TeV, A = 5 TeV, and no form factor (A = oo). 
The plot illustrates the great potential of the LHC to improve on existing bounds (even more so, when 
the machine will run at yet higher energy). 

For the production of four charged leptons, a similar NLO study (again including all off-shell, 
spin-correlation, virtual-photon-exchange, and interference effects) can be performed with the aMCONLO 
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Fig. 15: The total cross-section for Z, W + , and it, at the LHC (7 TeV) for NNPDF2.1 with a s (M z ) = 0.119 
(NLO and NNLO), a s (M z ) = 0.120 (NLO) and a s (M z ) = 0.114, 0.117 (NNLO), MSTW08 with a s (M z ) = 
0.1202 (NLO) and a s (M z ) = 0.1171 (NNLO), and ABKM09 with a s (M z ) = 0.1135 ± 0.0014 (NNLO). Un- 
certainties shown correspond to lcr. The band corresponds to the combination of CMS and ATLAS measurements. 



Figure taken from ref . Ml 241 . 



generator H 15311 . The parton showering can be done either with HERWIG and (for the first time in this 
context) with Pythia 6. The 0(a 2 s ) contribution of the gg channel is also included directly here, as 
obtained from MadLoop fl46l . In the study of ref. 111531 several key distributions together with the cor- 
responding theoretical uncertainties are presented. A further theoretical improvement is that scale and 
PDF uncertainties are computed at essentially no extra CPU-time using re-weighting techniques. 

6 Parton distribution functions and a s 

PDFs are the second thing one can not live without, if one works on LHC physics. Huge effort is devoted 
today in understanding differences and improving the theoretical and statistical treatment of PDFs. This 
activity is reflected in new PDFs sets being released by various groups II154II . The main focus of all 
groups is now directed towards NNLO PDFs sets, an improvement in the treatment of heavy quarks, an 
introduction of flexible parametrization, a more dynamic tolerance, and, of course, towards the inclusion 
of more data in the fits. Discussions are ongoing that try to clarify whether discrepancies between 
different PDFs are due to the inclusion of different data sets. For instance, there is no full consensus on 
what impact of the Tevatron jet data has on gluon distributions at the LHC. 

Fig. [TBI shows the uncertainty on three LHC benchmark processes (Z, W + , and it from the left 
to the right) coming from using different PDFs or a different value of a s , at NLO and at NNLO. Dif- 
ferences are due to the inclusion of different data in the fits, due to a different methodology (e.g. the 
parametrization), due to a different treatment of heavy quarks, and due to a different default value of the 
coupling constant. In particular, it is remarkable how much benchmark processes depend on the value 
of a s . The preliminary 2011 average value of a s is a s = 0.1183 ± 0.0010 111551 . It is interesting to 
note that the value barely changed compared to the 2009 number (a s = 0.1184 ± 0.0007) [156], but 
that the uncertainty on it increased. This is due to the inclusion of new data in the fits which tend to 
move the average value in opposite directions. An open issue today, in the combination of the various 
measurements to produce a world average for a s , is the treatment of outliers that have very small errors. 
This is the case for the extraction of a s from thrust computed at N 3 LL including power corrections using 
SCET 111571 . for the number obtained from T-decays in 111581 . and for the hadronic event shapes in e + e~ 
collisions at OPAL using NNLO+NLLA theoretical predictions 11591 . just to quote the most important 
cases. 

New processes added to the world average since 2009 include inclusive jets at the Tevatron [160], 
the e + e~ 3-jet rate which is know to NNLO H1611 . and e + e~ — > 5 jets, which is now known at 
NLO 1116211 . The use of the Tevatron jet data in this context is particularly interesting: while the error 
on a s from this extraction is not particularly small, this measurement and the sensitivity of benchmark 
processes on the value of a s shown in Fig. [15] raises the question of whether it is possible to make com- 
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Fig. 16: Jet mass in W + jets for a jet with substructure compatible with WZ(~> bb) or WH (— > bb) events. Figure 
taken from ref. 11 1 661 . 

petitive measurements of a s at the LHC. The extraction of the value of the coupling constant at hadron 
colliders must take into account that PDFs themselves do depend on a s . A viable possibility then is to 
consider appropriate ratios (e.g. Wj {Z + (n + 1) jets) / (Z/ (W + n jets))). 

7 Jet algorithms 

Jet algorithms are the third thing you can not live without, if you work on LHC physics. For a long 
time, infrared (IR) unsafe algorithms were used at the Tevatron, with several "patches" introduced to 
minimize the effect of the IR-unsafety. At the LHC, both ATLAS and CMS have adopted as default 
the anti-/c t algorithm [163]. Given that this algorithm was proposed only three years ago, it shows 
how flexible experimentalists are today in adopting new, successful ideasjj Using this algorithm both 
collaborations have already explored scales up to 4 TeV and could place constraints on various BSM 
models, in particular those models that would give rise to a resonance in the Mjj distribution (such as 
massive coloured bosons, black-holes, . . . ). 

Other IR-safe algorithms like the Cambridge-Aachen or SISCone are in use as well. These are 
particularly useful for studies which exploit the fact that when a massive boosted object decays, it gives 
rise to a "fat jet" with a non-trivial jet-substructure. Looking at the internal structure of these jets using 
jet-grooming techniques like filtering, pruning or trimming has a huge potential for making discoveries 
"easier" I1164[ ll651. These techniques have a big gain in sensitivity over traditional methods, but one 
might lose many events when imposing strict kinematical cuts and requiring a boosted regime. The po- 
tential of these studies has been demonstrated in several examples II 164111 651 . However sophisticated jet 
studies are still a young field, and as of now there are no precise rules on how to make discoveries easier. 
What is impressive, is that even these very new techniques are already being used at the LHC. Fig. [TBI 
shows for instance the single hadronic jet mass in W+ jet events in a boosted regime, an observable 
relevant for WH(— > bb). In Fig. [16] the Z peak coming from WZ{-^ bb) is evident, these very first 
results seem therefore very promising in finding a possible peak to due Higgs production. 

3 A minor downside to this is that ATLAS and CMS use a different radius - the choices for ATLAS are 0.4 and 0.6, while 
for CMS they are 0.5 and 0.7. 
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8 Conclusions 



The physics programme at the LHC is very rich: it spans from most precise measurements, e.g. in 
the case of Drell-Yan, to searches with highest reach for new physics either direct searches, or indirect 
ones (e.g. through the potential presence of anomalous couplings). This experimental programme at 
the LHC is supplemented by robust theoretical predictions that include NLO QCD corrections, mixed 
NLO-QCD+EW corrections, NNLO, and resummation of logarithmically enhanced contributions. Fur- 
thermore, different calculations are merged using clever matching procedures that catch the best features 
of different calculations. From the theoretical community, there is a clear and successful effort to pro- 
duce predictions and public codes that have the flexibility required for today's sophisticated experimental 
analysis (e.g. including parton shower, decays with spin correlations, massive quark effects, etc.). 

Impressive results have already come out of the LHC, but this is certainly only the tip of the 
iceberg. After just one year of running at the LHC some measurements, e.g. in processes involving W/Z 
production, start being dominated by theoretical and parton density errors. It will therefore be a real 
challenge for theorists to keep up with the high experimental precision. 
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